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Light nuclei can be produced in the central reaction zone via coalescence in rel- 
ativistic heavy ion collisions. E864 at BNL has measured the production of ten 
stable light nuclei with nuclear number of A = 1 to A = 7 at rapidity y ~ 1.9 and 
Pt/A < 300MeV/c. Data were taken with a Au beam of momentum of 11.5 A 
GeV/c on a Pb or Pt target with different experimental settings. The invariant 
yields show a striking exponential dependence on nuclear number over ten orders of 
magnitudes with a penalty factor of about 50 per additional nucleon. This penalty 
factor is used to estimate the strange quark matter (strangelet) production in the 
baryon rich and strangeness enhanced environment. The measurements of the pro- 
duction of antiproton, antideuteron, hypernuclei (^HAH), and strongly unstable 
nuclear states ( 4 H, 5 Li, 5 Li*, 5 He) are presented as well. A model of local ther- 
mal equilibrium with radial flow at the kinetic freeze-out with a temperature of 
T = 112 ± lOMeV, chemical potential of ti S = 503 ± 20MeV and flow velocity 
of about y/v]_ ~ 0.5c seems to be able to describe the data in the gross scale 
with the exceptions of the production of antihypcrons and hypernuclei. The large 
antihyperon production and the extra penalty for hypernuclei production are quite 
surprising. 



1. Introduction 

Relativistic heavy ion collisions may create high energy density and high 
baryon density in the reaction zone. Light nuclei can be produced by the 
recombination of created or stopped nucleonsEP. This recombination process 
is called coalescence. Coalescence of nuclear clusters can be characterized by 
the penalty factor for a nucleon added to the nuclear cluster. This idea can 
be extended to hypernuclei to calculate the strangeness penalty factor. Both 
baryon and strangeness penalty factor are useful to estimate the production 
rate of strange quark matter in terms of coalescence or thermal production. 
Since the probability of coalescence of a particular nuclear system (d, 3 ffe, 
etc.) depends on the properties of the hadronic system formed as a result of 
the collision, the study of the coalescence process is useful in elucidating those 
properties. For example, in a coalescence model, the coalescence probability 
depends on the temperature, baryon chemical potential (essentially the baryon 
density), and the size of the system, as well as the statistical weight (degen- 
eracy) of the coalesced nucleus. The data reported shows evidence that the 
probability may also depend on the binding energy of the coalesced nucleus at 
the level of about 6MeV. Systematic study of the production of light nuclei is 



1 



limited by their low production rates in relativistic heavy ion collisions. E864 
is the only experiment which is able to measure the production of charged nu- 
clei with A > 4 produced in the central reaction zone. However, the nature of 
the coalescence process tells us that the higher the baryon number, the more 
sensitive the production rate is to the system's configuration. For example, the 
inverse slope of the transverse momentum distribution (effective temperature) 
of "heavy" light nuclei is very sensitive to the density profile and radial flow 
profile. Production of antinuclci will provide additional information about the 
equation of state of the system and the state of chemical or kinetic equilibrium. 

In this report, the production of p, n, d, 3 iJe, t, 4 He, 6 He, 6 Li, 7 Li and 
7 Be around rapidity y ~ 1 .9 and transverse momentum of p^j A<~ SOOALeV/c 
in 10% most central Au+Pt(Pb) collisions measured by E864& at BNL is 
presented together with the production of p, J, 5 He, 5 Li, 4 iJ, \H and the 
production limits on 5 Li* and \H. These measurements have significant im- 
pact on the strange quark matterEf search by several experiments in relativistic 
heavy ion collisions. They also provide information about the thermal equi- 
librium of the system and the detailed process of coalescence. The richness of 
the data in rapidity and centralityo will not be discussed here. 

2. E864 Apparatus 

E864B is an open-geometry and high-data rate apparatus designed to 
search for new forms of matter. It has large rapidity and transverse momen- 
tum coverage. The tracking system consists of 2 dipole magnets for rigidity 
measurement and particle selection, 3 hodoscopes and 2 straw tube chambers 
measuring velocity, charge and position. The mass resolution is about 3% in the 
region of interest in high field settings. There is an additional mass measurement 
from hadronic calorimeter, which provides good energy and TOF information. 
These signals are used for a high mass level two trigger as well. Mom-detailed 
descriptions of the apparatus can be found in other publications&Q LrtJ . 

First, we can measure single-particle spectra of particles, such as 6 He, 
7 Li, 7 Be, n, p, and . Because E864 is open-geometry configuration, we 
have multiple particles in one event. We can also do two-particle correlation 
to reconstruct the invariant mass of their decay parent. We use event-mixing 
to subtract the background. Since we do not identify the decay vertex, strong 
decays and weak decays appear the same technically. Last but not the least, 
from the large sample of events, we can look for rare products (d) and search 
for new forms of matter such as strange quark matter. 

3. Results and Discussions 

In this section, we will present the results of the single-particle spectra 
and apply a simplified coalescence model assuming local thermal equilibrium 
to describe the data in the gross scale. Systematic errors of the temperature 
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Figure 1: Eighteen hadrons measured by E864. Stable light nuclei ^ are selected from 
y = 1.9 and p T /A < 300MeV/c. p is at y = 1.9 and p T ~ oW- d is at y = 1.9 
and < pt < lGeV/c(prcliminary). Hypernuclei are from 1.6 < y < 2.6 with correc- 
tion of detector's acceptance(preliminary). Strongly unstable nuclei are from y = 1.9 and 
pr5;2GeVyc(preliminary). Lines are fits to the data assuming local thermal equilibrium 
(work in progress). See text for details. 



and chemical potential extracted from this model-dependent simplificationlj 
will be discussed at the end. Radial flow is then studied using the temperature 
obtained and the inverse slope parameters from the transverse spectra of the 
K + , p, n, d, 3 He and 4 He (preliminary). Then we will show the results on 
the production of the hypernuclei and address the strangeness penalty factor. 
At the end, we will present the limits on strangelet search and discuss the 
significance of the limits in the context of the coalescence model. Particularly, 
we will discuss the limits on H-d (bound state of H-dibaryon and deuteron) 
search. 

Figure ^ shows the invariant yields of 18 hadrons we have measured in 
a selected rapidity (y) and transverse momentum (pt) bin. We observe a 
striking exponential dependence of light nuclei production up to A = 7 over 
ten orders of magnitudes in the invariant yield. Here, we will briefly present 
the results on d production as an example. Invariant yields are calculated 
and presented in terms of cP N / ' (2"npTdpTdy) in units of GeV~ 2 c 2 per central 
collision. In case of d, we have a sampled collection of 1.4 x 10 10 events of 
10% most central interactions. A signal of 17.6 ± 7.5 in the rapidity range of 
1.8 < y < 2.2 and 4.6 ± 3.3 in 1.4 < y < 1.8 above background is observed. 
The yields are determined to be 3.5 ± 1.5(stat)±Q 9 5 (sys) x !Q- 8 GeV~ 2 c 2 and 
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3.7 ± 2.7(stat)^\'i l (sys) x 10~ s GeV~ 2 c 2 respectively!]. These measurements 
suggest that the coalescence parameters B2 and B2 are consistent with each 
otherO. In this report, we will assume that antinuclei and nuclei have the 
same temperature (T) and chemical potential (|/k_b|) at the kinetic freeze-out. 
We fit the production of the ten stable light nuclei, with an exponential as a 
function of nuclear (baryon) number. The best fitEt is 26/48 A_1 . If we take 
our p and correct it for the antihyperon feeddownEP at 98% confidence level 
(CL)in y = 1.9 nearpr — 0, we have the fugacitytJ (Eq.6.6) of A 2 =exp( 2 ^ g -)= 
p/p — 23.3/0.0038. The fugacity can also be obtained from d/d: A' 4 = d/d = 
0.56/(3.5 x 10~ 8 ). The difference between A and A' comes from the difference 
of the isospin abundance or n/p ratio in a thermal distribution. From our 
unique measurements! of n/p = 1.19 ±0.08, we have A' = 1.1 A. For simplicity, 
the average fugacity or chemical potential are taken as: Xb = V AA', and 
[iB = {pp + Mn)/2. In fact, the fugacity from p/p at 98% CL is within the 
error of that from d/d. If we consider the penalty factor coming from the 
possibility of the population of the state as in a thermal system and neglect all 
the other details(flow has built up before the coalescence process and there is 
not additional momentum redistribution of the cluster after-burner) , we have 
exp((mjv — fJ-B)/T) ~ 48 where mjy is the nucleon mass. In Figure we 
plot the three equations above in the phase space of chemical potential and 
temperature. From the crossing point of the lines from penalty factor and 
fugacity from p/p, we get T = 112 ± lOMeV and /i b = 503 ± 20MeV (eraors 
are estimates, work in progress). This is quite close to other measurementsE3ljl . 
In the non-relativistic limit which for the case of nuclei is a good approximation, 
the kinetic energy can be separated into thermal energy and flow energy as 
T e ff = T + ^mV 2 in two-dimension radial flow. Figure || is the effective 
temperature vs. mass up to A — 4 with T = WIMeV at m = 0. Data of 
p, n, d, 3 He are at y = 2.3 where E864 has its best coverage in transverse 
momentum. For the flow at y = 2.3, we have to scale T (y = 1.9) and inverse 
slope parameter of A = A(y = 2.1) down by 10% to 20% which was not done in 
the plot. Therefore, the flow velocity we get probably will be an overestimate. 

The fit is T eff = 112 + 117 x A(MeV) with a mean velocity ^Jv[ ~ 0.5. 

Now we can use the flow and rapidity distribution measured to estimate 
the systematic error introduced by the model-dependent assumptionEJ that 
the system is in local thermal equilibrium. We can compare the difference 
between the total yield and the yield in one specific rapidity andjnomentum. 
When the parametrization of the correction factor of A x from Eq.6.10 is 
used, the limits of penalty between 39 and 72 are obtained for 1 > x > ~ h- 
If we assume that the temperature distribution in rapidity is a Gaussian with 
<Tt = 1.1 from our neutron measurements^ and the effective temperature 
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Figure 2: Temperature vs. chemical po- 
tential. jUta square and triangle points 
arc from E3D'. T = 112 ± lOAfeV, ^t s = 
503 ± 20MeV with estimated uncertainty. 
Errors are not shown. See text for details 
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Figure 3: Radial flow from transverse mo- 
mentum spectra. T and inverse slope of 
4 ffe are at y = 1.9 and y = 2.1 respec- 
tively which have higher temperature and 
stronger flow. The curved line is to guide 
the eye. 



scales as T e // — 1-0 + A with Boltzmann distribution, then T e ff(y) oc (1.0 + 
A) x exp (— (y — y cm ) 2 /(2<t^)). We have characterized the yiekLin rapidity at 
Pt — by the concavity!]! b/a and extrapolated to other nucleitf with limited 
data points. We ftbtain b/a = 0.14, 0.92, 2.8, 2.4, 3.2, 7.7 for p,d, 3 He, t, 4 iJe, 
6 He respectivelyQ'Q] . The penalty factor becomes 25 when we integrate over 
the transverse momentum and rapidity of —0.6 < y — y cm < 0.6(projectile and 
target regions are not included to avoid fragmentations, dN/dy oc T e ff(T e ff + 
m ) dyprdpT \pt=o)' wmcn should probably be taken as the absolute lower limit 
because the integration assumes non-correlation between momentum and space 
in the extrapolation to high pt for high mass which is in contradict with the 
flow phenomenon; in addition, there is indication that the effective temperature 
does not scale linearly with mass and might saturate at high masses as shown 
in Figure ||. In conclusion, the penalty factor is about 48^^ with uncertainty 
dominated by model-dependent estimate. 

E864 has a sampled collection of about 9 x 10 10 events of 10% most central 
interactions for hypernuclei studies. We have analyzed 2/3 of the data and 
found a \H signal at 1.8-sigma level above the background after mixed-event 
background subtraction^ . The invariant yield is 2.6 ± 1.4 x 10~ 4 Gev~ 2 c 2 
over 1.6 < y < 2.6 and < pr < 2GeV/c. The strangeness penalty factor is 
A s = \H/ 3 He = 0.036 ± 0.019 with an extra penalty of Ig^- = 0.17 ± 0.09 
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For the \H, we do not observe any significant signal above background and we 
set a upper limit at 90% CL of 0.29 x l{)~ A Gev~ 2 c 2 . When the degeneracy of 
S(2 J+ 1) = 4 is taken into account, we have As = \H / A/ A He < 0.032 which is 
consistent with the case. This extra penalty may come from the quantum 
mechanical correctionO . From this, we can write the cluster production^ as 
157/50" 4_1 /30' s '. The penalty for antimatter coalescing is much larger (~ 10 6 ). 

We have searched for strangelets (small strange quark matter) with charge 
of 0,±1, ±2, ±3 with null result; limits are set between 10 -9 to 10 -8 per central 
collision over a large mass range. We can discuss the sensitivity in terms of a 
coalescence mechanism and write the sensitivity a as A + 0.87|S'|$7.1 with the 
newly measured |Ag|. For example, we set a limit of 0.9 x 10 -7 for long-lived 
H-d state which has been predicted to be produced at a level of — * 10~ 7 
by C. Dover or 157/50 3 /30 2 = 1.4 x 10~ 6 from our measurementsB if it exists. 

In summary, we have measured the production of 18 hadrons in relativistic 
heavy ion collisions at the AGS energies. Information of the colliding system, 
such as temperature, chemical potential, flow and cluster production rate can 
be extracted from these measurements. Z. Xu would like to thank Dr. U. 
Heinz and Dr. J. Stachel for valuable discussions. We gratefully acknowledge 
the excellent support of the AGS staff. This work was supported in part by 
grants from the U.S. Department of Energy's High Energy and Nuclear Physics 
Divisions, and the U.S. National Science Foundation. 
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